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COMPUTER  CODE 


The  computer  code  used  in  the  investigation,  which  is 
described  in  References  1  and  2,  uses  the  implicit/explicit 
method  of  MacCormack  (Reference  3)  and  a  cartesian  grid.  Several 
different  turbulence  models  are  used,  depending  on  the  applica¬ 
tion.  In  some  of  the  studies,  where  only  averaged  quantities  are 
of  interest,  a  k-e  model  is  used  whereas  for  the  Very  Large  Eddy 
Simulations  (VLES)  a  vorticity  model  is  used  for  the  subgrid 
scales.  The  VLES  method  resolves  the  largest  scales  in  the  flow 
and  is  therefore  a  time  dependent  calculation.  The  method  is 
ideally  suited  to  the  present  range  of  problems  because  the  jet 
flows  are  dominated  by  the  streamwise  velocity  component  which 
has  a  peak  in  the  turbulent  kinetic  energy  spectrum  at  a  lower 
wave  number  than  the  other  components.  This  means  that  a  rela¬ 
tively  coarse  grid  will  resolve  a  reasonable  portion  of  the  large 
energy  containing  eddies. 

SPREADING  RATE  OF  THE  FOUNTAIN 


In  this  problem  two,  two-dimensional  wall  jets  colliding 
with  each  other  are  modeled.  At  the  jet  exits  a  suitable  repre¬ 
sentation  of  the  turbulence  is  used  and  this  is  validated  by 
comparing  computed  Reynolds  stresses  for  an  isolated  wall  jet 
with  experimental  data.  The  VLES  calculation  did  produce  the 
increased  spreading  rate  compared  with  an  averaged  calculation. 
The  mechanism  of  the  enhanced  spreading  rate  is  that  eddies  of 
different  kinetic  energy  collide  and  thus  the  high  energy  eddy 
will  push  through  the  opposing  flow  before  being  deflected  up¬ 
wards  in  the  fountain.  In  an  averaged  calculation  all  of  the 
colliding  eddies  have  the  same  energy  and  hence  do  not  push 
through  the  opposing  flow  but  are  deflected  upwards  about  the 
plane  of  symmetry.  This  work  is  reported  in  Reference  2. 
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The  accuracy  of  the  explanation  of  the  increased  spreading 
rate  given  above  is  enhanced  by  the  results  of  two  series  of 

calculations,  namely  the  introduction  of  heat  into  one  jet  as  a 
tracer  and  the  introduction  of  a  fence  at  the  center  of  the  col¬ 
lision  zone  of  the  two  jets.  The  VLES  technique  is  used  in  these 
studies. 

In  Figure  1  the  results  given  by  heating  the  right  wall  jet 
are  shown.  The  figure  shows  the  location  of  the  hot  gas,  at  some 
instant  and  it  can  be  seen  that,  the  fluid  from  the  right  hand 
jet  has  crossed  the  centerline.  Further  evidence  can  be  seen  in 
Figure  2  where  instantaneous  streaklines  are  shown.  The  streak¬ 
lines  in  both  Figure  2a  and  Figure  2b  are  at  a  location  cor¬ 
responding  to  the  wall  jet  thickness.  In  Figure  2a  the  streak¬ 
lines  from  the  left  hand  jet  are  shown  and  in  Figure  2b  those 

from  the  right  hand  jet  are  shown.  It  can  be  seen  that  the  fluid 
from  the  right  hand  jet  has  propagated  further  past  the  center- 
line  than  that  from  the  left  hand  jet,  corroborating  the  results 
of  the  heated  jet  case  shown  in  Figure  1. 

If  the  mechanism  for  the  increased  spreading  rate  is  that 
given  above  then  the  introduction  of  a  fence,  at  the  centerline, 
should  reduce  the  spreading  rate  since  fluid  from  the  jets  is 
inhibited  from  crossing  the  centerline  until  the  top  of  the 
fence.  Mach  number  contours  for  flows  with  and  without  the  fence 
are  shown  in  Figure  3  and  velocity  vectors  are  shown  in  Figure 
4.  It  can  be  seen  that  the  presence  of  the  fence  reduces  the 

spreading  rate  of  the  fountain  considerably. 

EFFECTS  OF  HEAT 

As  part  of  the  contract  the  effect  of  temperature  on  the 

flow  of  an  impinging  jet  was  examined  numerically.  In  these 
calculations  the  jet  was  at  a  temperature  of  1008  degrees  K  and 
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the  ambient  air  is  at  288  degrees  K.  The  jet  is  exhausting  into 
a  crossflow.  The  object  is  to  detmine  the  effect  of  temperature 
on  the  shape  of  the  ground  vortex  in  the  belief  that  this  is  a 
fairly  critical  parameter  in  such  flows.  The  calculations  are  a 
standard  version  of  a  high  Reynolds  number  k-e  turbulence  model 
for  a  variable  density  flow  described  by  Viegas  and  Rubesin  (Ref¬ 
erence  4);  wall  functions  are  used.  The  turbulent  kinetic 
energy,  k,  and  the  disappation  rate,  e,  are  given  by  the  follow¬ 
ing  equations.  The  turbulent  coefficient  of  viscosity  is  denoted 
by  uT.  The  molecular  viscosity  is  temperature  dependent  and 
satisfies  Sutherlands  law. 
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Two  runs  with  different  temperatures  but  with  similar  levels  of 
turbulent  kinetic  energy  at  the  jet  were  made  and  very  little 
effect  of  jet  temperature  was  observed.  It  should  be  noted  that 
variable  density  turbulence  models  have  not  been  tested  as  ade¬ 
quately  as  those  for  incompressible  flow  and  therefore  some 
reservation  should  be  attached  to  the  present  conclusion. 
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REYNOLDS  NUMBER  SCALING  OF  THE  SUCK  DOWN  EFFECT 

It  has  been  observed  that  small  scale  experiments  do  not 
produce  the  suck  down  effect  of  a  full  scale  test.  The  reasons 
for  the  discrepancy  are  not  known  although  there  are  several 
candidate  explanations,  such  as  Reynolds  number  effects  or  ini¬ 
tial  turbulence  levels  on  the  jet.  The  present  investigation  is 
concerned  with  possible  Reynolds  number  effects. 


The  difference  between  flows  at  different  Reynolds  number  is 
the  range  of  turbulence  scales  and  the  VLES  method  is  reasonably 
well  suited  to  modeling  a  range  of  scales.  In  a  VLES  calculation 
the  subgrid  scale  eddies  are  modeled  by  a  simple  vorticity  tur¬ 
bulence  model  and  the  eddy  size  that  is  resolved  is  related  to 
the  grid  size.  If  a  comparison  of  flows  at  two  different 
Reynolds  numbers  is  to  be  made  the  effect  of  the  range  of  scales 
is  captured  only  if  the  size  of  the  eddies  resolved  has  the  same 
relation  to  the  dissipation  length  scale  in  both  cases. 
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filtered  and  Reynolds  averaged  momentum 
calculation  is 
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equation  used  in 
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where  an  overbar  denotes  a  Reynolds  averaged  quantity  and  a  tilde 
denotes  a  quantity  that  is  resolved  by  the  grid.  The  analysis 
leading  to  this  equation  is  similar  to  that  given  by  Moin  & 
Rogallo  (Reference  5).  The  flow  quantities  are  scaled  to  jet 
exit  or  ambient  air  values  and  the  boundary  conditions  for  the 
impinging  jet  problem  are  invariant  with  Reynolds  number.  If  it 
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is  assumed  that  the  effect  of  molecular  viscosity  is  negligible 
then  the  only  Reynolds  number  effect  is  in  the  subgrid  scale  term 

C  A  2 
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where  Cg  is  a  constant  and  A  is  the  characteristic  size  of  the 


resolvable  eddies. 


If  the  size  of  the  resolved  scale  to  the 


Kolmogorov  scale  is  kept  constant  then  calculations  for  two  dif¬ 
ferent  Reynolds  numbers  will  model  the  difference  in  the  range  of 
scales . 
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The  Kolmogorov  length  scale,  n,  is  given  as 


n  =  ( 


where  U  is  a  characteristic  velocity  of  the  turbulence,  £  is  a 
characteristic  length  of  the  energy  containing  eddies  and  v  is 
the  kinematic  viscosity.  The  characteristic  velocity  can  be 
taken  to  be  the  velocity  of  the  jet  and  the  characteristic  length 
of  the  energy  containing  eddies  can  be  taken  to  be  the  VLES 
resolvable  grid  size  A,  since  A  will  be  at  a  particular  point  on 
the  turbulent  kinetic  energy  spectrum.  Thus,  in  order  to  repre¬ 
sent  the  range  of  scales,  the  ratio 


A  =  /UA) 

n  v 


must  be  a  constant.  Hence 


^  Li  ,  , 

—  =  constant 
v 

If  the  Reynolds  number  is  changed  by  changing  only  v  then  the 
grid  size  is  inversely  proportional  to  the  Reynolds  number.  Thus 
if  the  Reynolds  number  is  doubled  the  grid  size  is  halved. 
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Calculations  for  Reynolds  numbers  of  20,000,  40,000  and  60,000 
were  performed  for  the  configuration  shown  in  Figure  5,  namely 
two  wall  jets  colliding. 
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Since  the  grid  size,  A,  in  Eguation  1  is  multiplied  by  the 
constant  Cg  the  calculations  used  the  same  grid  which  makes  the 
truncation  error  uniform,  and  changes  Cg  in  inverse  proportion  to 
the  Reynolds  number. 

The  results  for  the  total  entrainment  are  given  below. 

Re:  20,000  40,000  60,000 

Entrainment:  0.345  0.419  0.370 

It  can  be  seen  that  there  is  an  effect  of  Reynolds  number  for 
Reynolds  numbers  of  this  order  of  magnitude. 

STUDIES  OF  FLOW  RESONANCE 

This  investigation  was  undertaken  to  determine  if  the  VLES 
technique  could  model  flow  resonance  such  as  that  studied  by  Ho 
and  Nossier  (Reference  6).  The  flows  considered  are  variations 
on  a  circular  jet  impinging  on  a  ground  plate.  The  height/jet 
diameter  ratio  is  3.1  and  the  nominal  jet  Mach  number  is  0.9. 
The  jet  exits  out  of  a  top  plate.  The  jet  has  a  constant  stagna¬ 
tion  pressure. 

The  first  flow  is  with  the  ground  plate  removed,  in  other 
words,  a  free  jet.  The  pressure  near  the  jet  centerline  approxi¬ 
mately  3  jet  diameters  down  from  the  jet  exit  is  shown  in  Figure 
6  as  a  function  of  time.  It  can  be  seen  that  the  pressure  does 
oscillate  periodically  with  time.  In  order  to  determine  if  there 
is  a  characteristic  Strouhal  number  the  fluctuating  pressure  was 


Fourier  transformed  to  give  p  and  the  power  spectrum  plotted 
against  the  Strouhal  number  (St)  ,  fD/Ujet,  where  f  is  the  fre¬ 
quency  and  D  is  the  jet  diameter.  In  Figure  (7)  the  power  spec¬ 
trum,  G,  of  the  pressure  fluctuations,  given  by 

...2 

G  -  |p| 

It  can  be  seen  that  there  is  a  peak  at  a  Strouhal  number  of  about 
0.45. 

The  second  flow  is  for  the  impinging  jet  and  in  Figure  (8) 
the  pressure  oscillations  close  to  the  ground  can  be  seen.  The 
power  spectrum  is  shown  in  Figure  (9)  and  does  not  have  a  spike 
at  the  Strouhal  number  of  0.45  but  the  general  behavior  is  not  as 
pronounced  as  in  the  free  jet. 

Finally,  an  attempt  to  model  more  closely  the  experiments  of 
Ho  and  Nossier  was  made  by  making  the  top  wall  mimic  the  presence 
of  the  nozzle  structure.  Results  are  shown  in  Figures  (10)  and 
(11).  It  can  be  seen  that  the  pressure  does  fluctuate  as  in  the 
previous  case  but  power  spectrum  is  considerably  different. 
However,  there  is  little  agreement  with  the  results  of  Ho  and 
Nossier . 

It  can  be  concluded  from  these  results  that  the  VLES  method 
can  represent  resonant  flows.  However,  the  quality  of  the 
results  need  to  be  improved  before  the  method  can  be  considered  a 
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Figure  2b.-  Streaklines  for  the  Right  Hand  Jet. 
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Figure  3a.-  Mach  Number  Contours  for  Jet  without  a  Fence 


A». 


T*WS' 


i  x  t  »  i  t  f  r  r  r  r  r 


■\  s  s  s  ». HHtflfflHli i > /•  r  /  /  r  / 


S  V  V  V  ». 

v ».  >>  v  *. — 


««c  /  /  /  /  / 


v  v  k  ^ w y^/frriirriT^iky *t ^ ^  /■  / 

V  V  *-  v.  *V  >v  /'/X^// 1  ^  A  A  A-  A  A  A  A 

^  ►  r.  V  /  ^  V1  <>•  ^ 


Sketch  of  Jet  Configuration 
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Figure  9.-  Variation  of  Power  Spectrum  with  Strouhal  Number  (Impinging  Jet) 


Figure  10.-  Variation  of  Pressure  with  Time  (Modified  Impinging  Jet). 


RESEARCH  OBJECTIVES 


The  objective  of  the  present  work  is  to  investigate  certain 
complex  fluid  interactions  for  impinging  jets.  These  interac¬ 
tions  may  be  unsteady.  The  investigation  uses  the  techniques  of 
Very  Large  Eddy  Simulation  (VLES)  of  the  flow  in  which  the  large 
eddies  in  the  flow  are  resolved  and  the  smaller  eddies  modeled. 
Although  computer  restrictions  do  not  allow  sufficient  grid 
points  for  a  quantitative  solution  of  the  problem  a  qualitative 
solution,  which  gives  insight  into  the  physical  aspects  is  pos¬ 
sible.  The  problems  that  have  been  studied  are  noted  below. 

a )  Spreading  Rate  of  the  Fountain 

When  two  wall  jets  collide  an  upwash  fountain  appears  which 
is  an  important  phenomena  in  V/STOL  aircraft  design.  It  had  been 
found  in  previous  work  (Contract  F49620-82-C-0031 )  that  the 
Reynolds  averaged  Navier  Stokes  equation  would  give  solutions 
with  approximately  one  third  of  the  spreading  rate  observed  in 
experiment.  Even  Reynolds  Stress  Transport  Models  (RSTM),  the 
most  advanced  turbulence  model  available,  did  not  resolve  this 
discrepancy  between  calculation  and  experiment.  Part  of  the 
present  work  is  to  use  the  VLES  technique  to  investigate  this 
problem. 

b)  Effects  of  Heat 

The  effects  of  high  temperature  in  the  jet  during  impinge¬ 
ment  are  investigated  to  determine  the  differences  between  hot 
and  cold  jet  flows. 


c )  Reynolds  Number  Scaling  of  The  Suck  Down  Effect 

It  had  been  observed  that  the  such  down  effect,  caused  by 
entrainment,  seemed  to  be  dependent  on  Reynolds  number.  Since 
the  Reynolds  averaged  equations  are  essentially  independent  of 
Reynolds  number  it  was  decided  to  examine  the  question  of 
Reynolds  number  dependence  using  the  VLES  technique.  The  differ¬ 
ence  between  a  high  and  low  Reynolds  number  is  the  range  of  tur¬ 
bulence  scales  and  this  can  be  modeled  by  the  VLES  method. 

d )  Flow  Resonance 

It  had  been  observed  in  experiments  in  impinging  jet  flows 
that  a  resonant  condition  could  appear  in  which  the  large  struc¬ 
tures  in  the  jet  would  lock  into  a  resonant  frequency.  It  was 
decided  to  attempt  to  model  this  phenomena  by  the  VLES  technique, 
based  on  the  assumption  that  the  largest  turbulent  structures  are 
responsible  for  the  resonance. 

ACCOMPLISHMENTS 

The  VLES  method  described  below  used  a  simple  vorticity 
model  for  the  subgrid  scales.  A  typical  turbulent  profile  was 
used  as  input,  with  its  time  average  approximately  equal  to  the 
experimental  values.  The  method  works  because  the  flows  under 
consideration  are  dominated  by  the  streamwise  velocity  compo¬ 
nents,  which  have  a  lower  wave  number  at  the  turbulent  kinetic 
energy  peak  than  the  other  components.  This  allows  the  use  of  a 
coarser  computational  grid  for  accuracy  that  might  otherwise  be 
the  case. 


a )  Spreading  Rate  of  the  Fountain 

In  this  problem  two,  two-dimensional  wall  jets  colliding 
with  each  other  are  modeled.  At  the  jet  exits  a  suitable  repre¬ 
sentation  of  the  turbulence  is  used  and  this  is  validated  by 
comparing  the  computed  Reynolds  stresses  for  an  isolated  wall  jet 
with  experimental  data.  The  VLES  calculation  did  produce  the 
increased  spreading  rate  compared  with  an  averaged  calculation. 
The  mechanism  of  the  enhanced  spreading  rate  is  that  eddies  of 
different  kinetic  energy  collide  and  thus  the  high  energy  eddy 
will  push  through  the  opposing  flow  before  being  deflected  up¬ 
wards  in  the  fountain.  In  an  averaged  calculation  all  of  the 
colliding  eddies  have  the  same  energy  and  hence  do  not  push 
through  the  opposing  flow  but  are  deflected  upwards  about  the 
plane  of  symmetry.  This  work  is  reported  in  Reference  1. 

In  addition  to  the  above  calculation  the  code  was  run  with  a 
fence  at  the  collision  point.  Since  the  fence  will  stop  the  high 
energy  eddies  at  the  plane  of  symmetry  it  is  to  be  expected  that 
the  spreading  rate  would  decrease.  This  conjecture  was  supported 
by  the  computational  results. 

As  a  final  case  the  code  was  run  again  with  one  jet  heated 
to  act  as  a  tracer.  This  calculation  supported  the  explanation 
of  the  enhanced  spreading  rate.  In  addition  it  indicated  that 
the  heat  transfer  could  not  be  modeled  by  gradient  diffusion. 

b )  Effects  of  Heat 

In  these  calculations  a  hot  jet  (1008  deg.  K)  exhausting 
into  an  atmosphere  at  288  deg.  K  with  a  crossflow  is  modeled. 
The  main  object  was  to  determine  the  effect  of  temperature  on  the 
shape  of  the  ground  vortex.  The  calculation  used  the  Reynolds 
averaged  equations  with  a  variable  density  k-e  turbulence  model 
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as  a  base.  The  viscosity  was  temperature  dependent  using 
Sutherlands  law.  Two  runs  with  different  temperatures  were  made, 
the  second  calculation  being  for  a  jet  at  ambient  temperature. 
The  difference  in  most  of  the  flow  parameters  was  of  the  order  of 
5%.  Hence  the  effects  of  temperature  do  not  appear  to  be  large. 

c)  Reynolds  Number  Scaling  of  the  Suck  Down  Effects 

In  these  calculations  the  VLES  code  for  two  colliding  wall 
jets  was  run  at  three  different  Reynolds  numbers,  20,000,  40,000, 
60,000.  The  grid  was  chosen  so  that  the  filter  size  in  each 
calculation  was  the  same  ratio  to  the  Kolmogorov  scale.  This 
does  give  a  Reynolds  number  effect  proportional  to  1/Re2.  In 
these  calculations  the  mass  flux  across  the  boundaries  excluding 
the  wall  jets  is  determined.  It  is  found  that  the  entrainment 
for  Re  =  2  x  104  is  approximately  25%  lower  than  that  for  Re  =  4 
x  104  while  at  Re  =  6  x  104  the  entrainment  is  approximately  8% 
lower  than  that  for  Re  =  4  x  104.  The  latter  change  may  be  due 
to  numerical  error.  However,  it  does  appear  that  there  is  a 
Reynolds  number  scaling  of  the  Reynolds  number  of  the  order  of 
104.  It  is  not  clear  whether  this  is  the  dominant  effect  in  the 
discrepancy  between  large  and  small  scale  experiments. 

d )  Resonance 

In  these  VLES  calculations  the  jet  has  a  Mach  number  of  0.9 
and  a  height/diameter  of  3.  Several  attempts  were  made  to  induce 
resonance  but  no  definite  conclusion  was  reached.  It  did  appear 
as  if  a  resonance  phenomena  was  happening  in  bursts  but  it  is  not 
clear  if  this  is  due  to  numerical  or  physical  aspects  of  the 
calculation.  A  calculation  of  a  free  jet  did  develop  a  periodic 
structure  of  ring  vortices  which  indicates  that  the  method  can 
model  large  periodic  structures.  An  interesting  aspect  that  did 
arise  in  the  calculations  was  the  appearance  of  shock  waves  near 
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the  impingement  zone  due  to  the  turbulence.  These  shock  waves 
appeared  either  at  the  edge  or  at  the  top  of  the  impingement 
zone . 
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